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Abstract. Wc report the angular diameter measurements of seven classical Cepheids, XSgr, 7?Aql, WSgr, ("Gem, 
(5 Dor, Y Oph and £ Car that we have obtained with the VINCI instrument, installed at ESO's VLT Interferometer 
■ (VLTI). We also present reprocessed archive data obtained with the FLUOR/IOTA instrument on £ Gem, in 

order to improve the phase coverage of our observations. We obtain average limb darkened angular diameter 
values of 0ZrJ[XSgr] = 1.471 ± 0.033 mas, Ko[r] Aq\] = 1.839 ± 0.028 mas, 0ZrJ[WSgr] = 1.312 ± 0.029 mas, 
6>LD[/3Dor] = 1.891±0.024mas, 6> L D[CGem] = 1.747±0.061 mas, L D[YOph] = 1.437±0.040 mas, and LD [ICar] = 
2.988 ± 0.012 mas. For four of these stars, r\ Aql, WSgr, /3Dor, and ^Car, we detect the pulsational variation of 
their angular diameter. This enables us to compute directly their distances, using a modified version of the Baade- 
Wesselink method: d [77 Aql] = 276±|| pc, d[WSgr] = 379±?^ pc, d[/3Dor] = 345+^ 5 pc, d[£Cax] = 603^ 4 pc. 
The stated error bars are statistical in nature. Applying a hybrid method, that makes use of the Gieren et 
al. (1998) Period- Radius relation to estimate the linear diameters, we obtain the following distances (statistical 
and systematic error bars are mentioned): d[X Sgr] = 324±7±17pc, d[r? Aql] = 264±4±14 pc, d[W Sgr] = 386±9± 
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21 pc, d[/3Dor] = 326±4±19pc, d[CGem] = 360±13±22 pc, d[Y Oph] = 648±17±47pc, d[£ Car] = 542±2±49pc. 
Key words. Techniques: interferometric, Stars: variables: Cepheids, Stars: oscillations 
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1. Introduction authors (see e.g. Benedict et al. 20U^. that the LMC dis- 

tance is currently the weak link in the extragalactic dis- 

bor almost a century, Cepheids have occupied a cen- , , , , , A . , , , . 

J \ r . . 1 • • 1 1 tance scale ladder. Another avenue is to determine the 

tral role in distance determinations. This is thanks to r ,i n T ^ ^ t~\ i i -j 

zero point ol the P-L relation with Galactic Cepheids, us- 

the existence ol the Period-Luminosity (P-L) relation. f . , n , ^ , . , 

; ' ' mg lor instance parallax measurements, Cepheids m cms- 

M = a l0gP + b ; w ^ lch / el j ates the logarithm of the van- tergj Qr t h thg Baade _ Wesselink (BW) method . We 

ability period ol a Cepheid to its absolute mean magm- .1 ■ 1 1 m TT j ttt\ 

J , , . propose m this paper and its sequels (Papers 11 and III) 

tude. ihesc stars became even more important since the , . ,. r ,-, r , , . , n ^ n T , 

TT , , . * to improve the calibration ot the Cepheid P-R, P-L and 

Bubble Space lelescope Key Project on the cxtragalac- r , . , , , , ,. ,-, , , . ,. 

. r , J ,11. , surtace brightness-color relations through a combination 

tic distance scale (Inccdman et al. 2001) has totally relied c j - , r u n. • 1. 

; ~~ — ' ol spectroscopic and mterierometric observations ol bright 

on Cepheids for the calibration of distance indicators to , ,. ^ , . , 

, , „ Galactic Cepheids. 

reach cosmologically significant distances. In other words, . 

■t 4-1, ru 4.- c ii. n i, • j r> t 1 t,- ■ In the particular case 01 the P-L relation, the slope 

it the calibration ot the Cepheid P-L relation is wrong, . / ^ t „.--., , ^ \ ■ , , 

, . a is well known Irom Magellanic Cloud Cepheids (e.g. 

the whole extragalactic distance scale is wrong. TT . . _l v 

. T-i .1 T-i t 1 Udalski et al. 1999), though Lanorx et al. ( 1999) have sug- 

lhere are various ways to calibrate the P-L rela- , : . ~ , , , * . , , 

, , T Trtrn t/ n • gested that a Malmquist etlect (population mcomplete- 

tion. I he avenue chosen by the Hbl Key-Project was . , , , . , . , ^ , , , , ,. 

,. , .-1 -[- , , . ™ , ness) could bias this value. On the other hand, the can- 
to assume a distance to the Large Magellanic Cloud . . . . ■,,/,, , , • 1 \ 
/t urn 1 j " ■ ■ c ia j- 1 bration ot the zero-pomt b the hypothetic absolute mag- 
(LMC), thereby adopting a zero point ot the distance . , ^ , ■ ,\ 

. ' . , , . JW ,n .. . ,. mtude ot a 1-day period Cepheid) requires measurement 
scale, .breedman et al. (120011) present an extensive discus- .. / . . ,>,,., .„ , . , 
r .... T ,vZ ' ' , ... ot the distance to a number ot nearby Cepheids with high 
sion ot all available LMC distances, and note, with other . r . . 
precision, .tor this purpose, mterterometry enables a new 
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star's temperature, as we have directly access to its an- 
gular diameter fDavi sTH)79l Sasselov & Karovska 1994). 
Using this method, we derive directly the distances to the 
four nearby Cepheids 77 Aql, WSgr, ,3 Dor and £Car. For 
the remaining three objects of our sample, X Sgr, £ Gem 
and Y Oph, we apply a hybrid method to derive their dis- 
tances, based on published values of their linear diameters. 

After a short description of the VINCI/VLTI instru- 
ment (Sect.EJ, we describe the sample Cepheids that we 
selected (Sect. OH)- In Sect, ri] and 03 we report our new ob- 
servations as well as reprocessed measurements of £ Gem 
retrieved from the FLUOR/IOTA instrument archive. 
Sect.||3is dedicated to the computation of the correspond- 
ing angular diameter values, taking into account the limb 
darkening and the bandwidth smearing effects. In Sect.0 
and |H1 we investigate the application of the BW method 
to our data, and we derive the Cepheid distances. 

We will discuss the consequences of these results 
for the calibration of the Period- Radius (P-R), Period- 
Luminosity (P-L) and Barnes-Evans relations of the 
Cepheids in forthcoming papers (Paper II and III). 

2. Instrumental setup 

The European Southern Observatory's Very Large 
Telescope Interferometer (Glindemann et al. 2000) is in 
operation on Ccrro Paranal, in Northern Chile since 
March 2001. For the observations reported in this pa- 
per, the beams from two Test Siderostats (0.35 m aper- 
ture) or two Unit Telescopes (8 m aperture) were re- 
combined coherently in VINCI, the VLT INterferometer 
Commissioning Instrument (Kervella et al. 120001 12003a). 
We used a regular K band filter (A = 2.0 — 2.4 /im) that 
gives an effective observation wavelength of 2.18 (im for 
the effective temperature of typical Cepheids (see Sect. 16.41 
for details). Three VLTI baselines were used for this pro- 
gram: E0-G1, B3-M0 and UT1-UT3 respectively 66, 140 
and 102.5 m in ground length. Fig.dshows their positions 
on the VLTI platform. 

3. Selected sample of Cepheids 

Despite their brightness, Cepheids are located at large 
distances, and the Hipparcos satellite (Perryman et 
al. 119971) could only obtain a limited number of Cepheid 
distances with a relatively poor precision. If we exclude 
the peculiar first overtone Cepheid aUMi (Polaris), the 
closest Cepheid is JCep, located at approximately 250 pc 
(Mourard et al. 119971 Nordgren et al. l2000[> . As described 
by Davis l|1979fl and Sasselov & Karovska l|1994fl . it is pos- 
sible to derive directly the distance to the Cepheids for 
which we can measure the amplitude of the angular di- 
ameter variation. Even for the nearby Cepheids, this re- 
quires an extremely high resolving power, as the largest 
Cepheid in the sky, £Cai, is only 0.003" in angular di- 
ameter. Long baseline interferometry is therefore the only 
technique that allows us to resolve these objects. As a re- 
mark, the medium to long period Cepheids (D 200 D Q ) 




Fig. 1. Layout of the three baselines used for the 
VINCI/VLTI Cepheids observations, UT1-UT3 (102.5m), 
E0-G1 (66m) and B3-M0 (140m). 



in the Large Magellanic Cloud (LMC) (d « 55kpc) are 
so small (9 ~ 30 /xas) that they would require a baseline 
of 20 km to be resolved in the K band (5 km in the vis- 
ible). However, such a measurement is highly desirable, 
as it would provide a precise geometrical distance to the 
LMC, a critical step in the extragalactic distance ladder. 

Mourard ( 1996) has highlighted the capabilities of the 
VLTI for the observation of nearby Cepheids, as it pro- 
vides long baselines (up to 202 m) and thus a high resolv- 
ing power. Though they are supergiant stars, the Cepheids 
are very small objects in terms of angular size. A conse- 
quence of this is that the limit on the number of interfer- 
omctrically resolvable Cepheids is not set by the size of 
the light collectors, but by the baseline length. From pho- 
tometry only, several hundred Cepheids can produce in- 
terferometric fringes using the VLTI Auxiliary Telescopes 
(1.8m in diameter). However, in order to measure accu- 
rately their size, one needs to resolve their disk to a suffi- 
cient level. Kervella <|2001afl has compiled a list of more 
than 30 Cepheids that can be measured from Paranal 
using the VINCI and AMBER (Petrov et al. 2000) in- 
struments. Considering the usual constraints in terms of 
sky coverage, limiting magnitude and accessible resolu- 
tion, we have selected seven bright Cepheids observable 
from Paranal Observatory (latitude A = — 24deg): XSgr, 
77 Aql, WSgr, /3Dor, (Gem, YOph and £ Car. The pe- 
riods of these stars cover a wide range, from 7 to 35.5 
days. This coverage is important to properly constrain 
the P-R and P-L relations. To estimate the feasibility 
of the observations, the angular diameters of these stars 
were deduced from the BW studies by Gieren et al. (1993). 
For C Gem and r\ Aql, previously published direct interfer- 
ometric measurements by Nordgren et al. (2000), Kervella 
et al. (|2001b|l and Lane et al. (P002 ) already demonstrated 
the feasibility of the observations. The relevant parame- 
ters of the seven Cepheids of our sample, taken from the 
literature, are listed in Tabled 
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Table 1. Relevant parameters of the observed sample of Cepheids, sorted by increasing period. 
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X Sgr 


rj Aql 


W Sgr 


Dor 


C Gem 


Y Oph 


i Car 




HD 161592 


HD 187929 


HD 164975 


HD 37350 


HD 52973 


HD 162714 


HD 84810 


a 

my 


A COI 

4.D81 


O ClA O 

3.942 


4. iUU 


3. (31 


3.928 


6.164 


6.1 1 1 


m K b 


2.56 


1.966 


2.82 


1.959 


2.11 


2.682 


1.091 


Sp. Type 


F5-G2II 


F6Ib-G4Ib 


F4-G2Ib 


F4-G4Ia-II 


F7Ib-G3Ib 


F8Ib-G3Ib 


F6Ib-K0Ib 


n (mas) c 


3.03 ± 0.94 


2.78 + 0.91 


1.57 + 0.93 


3.14 + 0.59 


2.79 + 0.81 


1.14 + 0.80 


2.16 ± 0.47 


Min T cff (K) 


5670 


5400 


5355 


5025 


5150 






Mean T cff (K) d 


6150 


5870 


5769 


5490 


5430 


5300 


5090 


Max T cff (K) 


6820 


6540 


6324 


6090 


5750 






Min log q 


1.86 


1.25 


1.72 


1.60 








Mean log<? c 


2.14 


1.49 


1.82 


1.83 


1.50 


1.50 


1.50 


Max \ogg 


2.43 


1.73 


2.02 


2.06 








[M/H] e 


0.04 


0.05 


-0.01 


-0.01 


0.04 


0.05 


0.30 


To (JD-2.452.10 6 ) f 


723.9488 


519.2477 


726.8098 


214.2153 


210.7407 


715.4809 


290.4158 


P (days) g 


7.013059 


7.176769 


7.594904 


9.842425 


10.150967 


17.126908 


35.551341 


Intensity profiles 11 
















ai 


+0.7594 


+0.8816 


+0.8002 


+0.7969 


+0.8713 


+0.8549 


+0.8500 




-0.4530 


-0.7418 


-0.5135 


-0.4596 


-0.6536 


-0.5602 


-0.4991 


03 


+0.0347 


+0.3984 


+0.1583 


+0.1341 


+0.3283 


+0.2565 


+0.2113 


0,4 


+0.0751 


-0.0778 


+0.0109 


+0.0082 


-0.0610 


-0.0437 


-0.0340 



a m v from Barnes et al. H1987H for XSgr, from Barnes et al. •: 1997: for rj Aql, from Moffett & Barnes J1984H for WSgr and 
C Gem, from Berdnikov & Turner l^rJTjTl for Dor and I Car, and from Coulson & Caldwell tTMEi for Y Oph. 

b rriK from Welch et al. fi'gfg4^ for XSgr, and WSgr, from Laney & Stobie 1199211 for 0~Dor, YOph, and IGax, from Ducati et 
al. ll2Wfl for C Gem, from Barnes et al. ftWl for t) Aql. 

c Parallaxes from the HlPPARCOS catalogue (Perryman et al. 119971 . 

d From Kiss & Szatmarv 113551 for (Gem and jyAql, Bersier et al. (TlWI for WSgr, and Pe HtT?75l for XSgr and /3Dor. 
6 From Andrievsky et al. lT20U2Tl . Cayrel de Strobel et al. (TTWI IStMl . and Pel ltT?78l . except for log 5 of YOph. 
f Reference epoch To values have been computed near the dates of the VINCI observations, from the values published by 
Szabados lU9B9aL 

g P values from Szabados ( 1989a I . The periods of r\ Aql, £ Gem and W Sgr are known to evolve. The values above correspond 
to the To chosen for these stars. 

Four-parameters intensity profiles from Claret (TOOO) in the K band, assuming a microturbulence velocity of 4km/s and the 
average values of T e ff and log g. 



4. Interferometric data processing 

4.1. Coherence factors 

We used a modified version (Kervella et al. l2003c[) of 
the standard VINCI data reduction pipeline, whose gen- 
eral principle is based on the original algorithm of the 
FLUOR instrument (Coude du Foresto et al. 119971 Coude 
du Foresto et al. H998a|) . The VINCI/VLTI commissioning 
data we used for this study are publicly available through 
the ESO Archive, and result from two proposals of our 
group, that were accepted for ESO Periods 70 and 71. 

The goal of the raw data processing is to extract the 
value of the modulated power contained in the interfero- 
metric fringes. This value is proportional to the squared 
visibility V 2 of the source on the observation baseline, 
which is in turn directly linked to the Fourier transform 
of the light distribution of the source through the Zernike- 
Van Cittert theorem. 

One of the key advantages of VINCI is to use single- 
mode fibers to filter out the perturbations induced by the 



turbulent atmosphere. The wavefront that is injected in 
the fibers is only the mode guided by the fiber (gaussian 
in shape, see Ruilierfl.999 or Coude du Forestc fT998bl for 
details) . The atmospherically corrupted part of the wave- 
front is not injected into the fibers and is lost into the 
cladding. Due to the temporal fluctuations of the turbu- 
lence, the injected flux changes considerably during an 
observation. However, VINCI derives two photometric sig- 
nals that can be used to subtract the intensity fluctua- 
tions from the interferometric fringes and normalize them 
continuously. The resulting calibrated interferograms are 
practically free of atmospheric corruption, except the pis- 
ton mode (differential longitudinal delay of the wavefront 
between the two apertures) that tends to smear the fringes 
and affect their visibility. Its effect is largely diminished 
by using a sufficiently high scanning frequency, as was the 
case for the VINCI observations. 

After the photometric calibration has been achieved, 
the two interferograms from the two interferometric out- 
puts of the VINCI beam combiner are subtracted to re- 
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move the residual photometric fluctuations. As the two 
fringe patterns are in perfect phase opposition, this sub- 
traction removes a large part of the correlated fluctu- 
ations and enhances the interferometric fringes. Instead 
of the classical Fourier analysis, we implemented a time- 
frequency analysis (Segransan et al. I1999f) based on the 
continuous wavelet transform (Farge 1992). In this ap- 
proach, the projection of the signal is not onto a sine wave 
(Fourier transform), but onto a function, i.e. the wavelet, 
that is localised in both time and frequency. We used as a 
basis the Morlet wavelet, a gaussian envelope multiplied 
by a sine wave. With the proper choice of the number of os- 
cillations inside the gaussian envelope, this wavelet closely 
matches a VINCI interfcrogram. It is therefore very effi- 
cient at localizing the signal in both time and frequency. 

The differential piston corrupts the amplitude and the 
shape of the fringe peak in the wavelet power spectrum. 
A selection based on the shape of fringe peak in the time- 
frequency domain is used to remove "pistonned" and false 
detection interferograms. Squared coherence factors /i 2 are 
then derived by integrating the wavelet power spectral 
density (PSD) of the interferograms at the position and 
frequency of the fringes. The residual photon and detector 
noise backgrounds are removed by making a least squares 
fit of the PSD at high and low frequency. 

4.2. Calibrators 

The calibration of the Cepheids' visibilities was achieved 
using well-known calibrator stars that have been selected 
in the Cohen et al. lfT?9 9) catalogue, with the exception of 
elnd. This dwarf star was measured separately (Segransan 
et al. l2003[) and used to calibrate one of the r\ Aql measure- 
ments. The angular diameters of 39EriA, HR4050 and 
HR 4546 (which belong to the Cohen et al. H999l catalogue) 
were also measured separately, as these stars appeared to 
give a slightly inconsistent value of the interferometric ef- 
ficiency. 

For 39 Eri A and HR 4546, the measured angular diam- 
eters we find are 6» UD = 1.74 ± 0.03 and 2.41 ± 0.04 mas, 
respectively. These measured values are only 2 a lower 
than the Cohen et al l|1999|l catalogue values of #ud = 
1.81 ± 0.02 and 2.53 ± 0.04 mas. A possible reason for this 
difference could be the presence of faint, main sequence 
companions in orbit around these two giant stars. The 
additional contribution of these objects would bias the di- 
ameter found by spectrophotometry towards larger values, 
an effect consistent with what we observe. For HR4050, we 
obtained 6>tjd = 5.18 ± 0.05 mas, only +1 a away from the 
catalogue value of #ud = 5.09 ±0.06 mas, . The character- 
istics of the selected calibrators are listed in Table |21 The 
limb-darkened disk (LD) angular diameters of these stars 
were converted into uniform disk values using linear coef- 
ficients taken from Claret et al. ©95). As demonstrated 
by Borde et al. I|2002[) . the star diameters in the Cohen 
et al. ( 1999 ) list have been measured very homogeneously 



to a relative precision of approximately 1% and agree well 
with other angular diameter estimation methods. 

The calibrators were observed soon before and after 
the Cepheids, in order to verify that the interferomet- 
ric efficiency (IE) has not changed significantly during 
the Cepheid observation itself. In some cases, and due to 
the technical nature of commissioning observations, part 
of the Cepheid observations could not be bracketed, but 
only immediately preceded or followed by a calibrator. 
However, the stability of the IE has proved to be gener- 
ally very good, and we do not expect any significant bias 
from these single-calibrator observations. Some observa- 
tions included several calibrators to allow a cross-check 
of of their angular sizes. The calibrators were chosen as 
close as possible in the sky to our target Cepheids, in or- 
der to be able to observe them with similar airmass. This 
selection has taken into account the constraints in terms 
of limiting magnitude and sky coverage imposed by the 
VLTI siderostats and delay lines. The IE was computed 
from the coherence factor measurements obtained on the 
calibrators, taking into account the bandwidth smearing 
effect (see Sect. 16.41) and a uniform disk angular diameter 
model. This calibration process yielded the final squared 
visibilities listed in Tables |3 to 



5. Data quality 

5.1. General remarks 

Due to the fact that we used two types of light collec- 
tors (siderostats and UTs) and several baselines (from 66 
to 140m in ground length), the intrinsic quality of our 
data is relatively heterogeneous. In this Section, we dis- 
cuss briefly the characteristics of our observations of each 
target. One particularity of our measurements is that they 
have all been obtained during the commissioning period of 
the VLTI, during which technical tasks were given higher 
priority. In particular, the long baseline B3-M0 was only 
available during a few months over the two years of opera- 
tions of the VLTI with VINCI. The UT1-UT3 observations 
were executed during two short commissioning runs and 
it was not possible to obtain more than one or two phases 
for the observed stars (/3Dor and CGem). However, the 
very large SNR values provided by the large aperture of 
the UTs, even without high-order adaptive optics, gave 
high-precision visibility measurements. 

The VINCI processing pipeline produces a number of 
outputs to the user for the data quality control, includ- 
ing in particular the average wavelet power spectral den- 
sity (WPSD) of the processed interferograms. This is an 
essential tool to verify that no bias is present in the cali- 
brated and normalized fringe power peak. Fig.[3shows the 
average WPSD of a series of 302 interferograms obtained 
on XSgr. No bias is present, and the residual background 
is very low. The power integration being done between 
2000 and 8000 cm -1 , the complete modulated power of 
the fringes is taken into account without bias. 
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Table 2. Relevant parameters of the calibrators. 



Name 






772. Y 




Sd TvDe 


- 1 - ctt V / 


lo 




7r (mas) a 


^Tnfmas 1 ) 13 




xPhe 


HD 


12524 


5.16 


1.52 


K5III 


3780 


1 


9 


8.76 ±0.64 


2.77 ±0.032 


2.69 ±0.031 


39 Eri A 


HD 


26846 


4.90 


2.25 


K3III 


4210 


2 


2 


15.80 ± 0.95 


1.79 ±0.031* 


1.74 ± 0.030* 


eRet 


HD 


27442 


4.44 


1.97 


K2IVa 


4460 


2 


3 


54.84 ± 0.50 


1.95 ± 0.049 


1.90 ± 0.048 


HR2533 


HD 


49968 


5.69 


2.10 


K5III 


3780 


1 


9 


6.36 ± 0.92 


1.93 ±0.020 


1.87 ±0.019 


HR2549 


HD 


50235 


5.00 


1.39 


K5III 


3780 


1 


9 


3.60 ± 0.56 


2.25 ± 0.036 


2.18 ±0.035 


7 2 Vol 


HD 


55865 


3.77 


1.52 


K0III 


4720 


2 


6 


23.02 ± 0.69 


2.50 ± 0.060 


2.44 ± 0.059 


6 Pup 


HD 


63697 


5.18 


2.62 


K3III 


4210 


2 


2 


12.87 ±0.71 


1.88 ± 0.039 


1.83 ±0.038 


HR3046 


HD 


63744 


4.70 


2.31 


K0III 


4720 


2 


6 


14.36 ± 0.48 


1.67 ±0.025 


1.63 ± 0.024 


HR4050 


HD 


89388 


3.38 


0.60 


K3IIa 


4335 


2 


3 


4.43 ± 0.49 


5.32 ± 0.050* 


5.18 ±0.048* 




TTPi 

riJJ 


onnno 

oyyyo 




Z.4U 


TV" 1 TTT 


a £ en 


2 


- 




lb. ZD ± U.DO 


1. 1 Z rfc U.UzU 


1 RQ _1_ A Al A 


HR4526 


HD 


102461 


5.44 


1.77 


K5III 


3780 


1 


9 


3.97 ±0.61 


3.03 ± 0.034 


2.94 ± 0.033 


HR4546 


HD 


102964 


4.47 


1.56 


K3III 


4210 


2 


2 


7.03 ± 0.72 


2.48 ± 0.036* 


2.41 ± 0.035* 


HR4831 


HD 


110458 


4.67 


2.28 


K0III 


4720 


2 


6 


17.31 ± 0.65 


1.70 ±0.018 


1.66 ±0.018 


xSco 


HD 


145897 


5.25 


1.60 


K3III 


4210 


2 


2 


7.43 ± 0.91 


2.10 ±0.023 


2.04 ± 0.022 


70Aql 


HD 


196321 


4.90 


1.21 


K5II 


3780 


1 


9 


1.48 ± 0.91 


3.27 ±0.037 


3.17 ±0.036 


7Aqr 


HD 


199345 


5.50 


2.00 


K5III 


3780 


1 


9 


5.42 ± 0.99 


2.14 ±0.024 


2.08 ±0.023 


eind 


HD 


209100 


4.69 


2.18 


K4.5V 


4580 


4 


5 


275.79 ± 0.69 


1.89 ±0.051* 


1.84 ±0.050* 


AGru 


HD 


209688 


4.48 


1.68 


K3III 


4210 


2 


2 


13.20 ± 0.78 


2.71 ± 0.030 


2.64 ±0.029 


HR8685 


HD 


216149 


5.41 


1.60 


M0III 


3660 


1 


4 


2.95 ± 0.69 


2.07 ±0.021 


2.01 ±0.020 



a Parallaxes from the HlPPARCOS catalogue (Perryman et al. 119971 . 

b Catalogue values from Cohen et al. |1999|, except for eind, HR4050 and 39 Eri A. 

c Linear limb darkening coefficients factors from Claret et al. |1995). 

* The angular diameters of eind, HR4050, HR4546 and 39 Eri A have been measured separately with VINCI. 




to its declination of 5 = —28 deg, X Sgr culminates almost 
at zenith over Paranal (-24 deg), and all the observations 
were obtained at very low airmasses. It is located on the 
sky near two other Cepheids of our sample, YOph and 
WSgr, and these three stars share the same calibrator, 
X Sco. The average signal to noise ratio (SNR) was typi- 
cally 2 to 5 on the photometric outputs of VINCI, and 4 to 
6 on the interferometric channels, for a constant fringe fre- 
quency of 242 Hz. A total of 4 977 interferograms were pro- 
cessed by the pipeline. The same remarks apply to W Sgr 
and Y Oph, as they have almost the same magnitude and 
similar angular diameters. The number of processed inter- 
ferograms for these two stars was 4 231 and 2 182, respec- 
tively, during 9 and 4 observing sessions. 



4000 6000 8000 

Wavenumber {cm A -1 ) 



Fig. 2. Average wavelets power spectral density of 302 
interferograms obtained on XSgr on JD=2452768.8462. 
No background or bias is present. The integration of 
the fringes modulated power is done between 2000 and 
8000 cm- 1 



5.2. X Sgr, WSgr and YOph 

X Sgr was observed 8 times on the B3-M0 baseline (140 m 
ground length), using exclusively the two 0.35m Test 
Siderostats (TS). The projected baseline length varied be- 
tween 118.4 and 139.7 m, and the observed squared visibil- 
ities were confined between V 2 = 56.9 and 71.1 %. Thanks 



5.3. r/Aql 



i] Aql was observed once on the E0-G1 baseline (66 m) and 
10 times on the B3-M0 baseline (140m ground length). 
The total number of processed interferograms is 5 584. 
The SNRs were typically 4 and 7 on the photometric and 
interferometric outputs at a fringe frequency of 242 to 
272 Hz. Due to its northern declination (6 = ±1 deg) and 
to the limits of the TS, it was not possible to observe 
r\ Aql for more than two hours per night, therefore limit- 
ing the number of interferograms and the precision of the 
measurements. 
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5.4. (3 Dor 

/3Dor is a difficult target for observation with the TS, 
as it is partially hidden behind the TS periscopes that 
are used to direct the light into the VLTI tunnels. This 
causes a partial vignetting of the beams and therefore a 
loss in SNR. The data from the TS are thus of interme- 
diate quality, considering the brightness of this star. It 
is located at a declination of —62 deg, relatively close to 
£ Car, and therefore these two stars share some calibrators. 
In addition to the 5 observations with the TS, four mea- 
surements were obtained during three commissioning runs 
on the UT1-UT3 baseline. A total of 8 129 interferograms 
were processed, of which 5 187 were acquired with the 8 m 
Unit Telescopes (96 minutes spread over four nights were 
spent on /3Dor using UT1 and UT3). 




100 120 140 160 
Projected baseline (m) 



5.5. ( Gem 

At a declination of +20 deg, £ Gem is not accessible to 
the TS due to a mechanical limitation. This is the reason 
why this star was observed only on two occasions with 
UT1 and UT3, for a total of 3 857 interferograms, obtained 
during 41 minutes on the target. The average on-source 
SNRs were typically 50 for the interferometric channels 
and 30 for the photometric signals, at a fringe frequency 
of 694 Hz. 

The data obtained using the FLUOR/IOTA instru- 
ment are described in Kervella et al. ( 2001b). They were 
reprocessed using the latest release of the FLUOR soft- 
ware that includes a better treatment of the photon shot 
noise than the 2001 version. As the baseline of IOTA is 
limited to 38 m, the visibility of the fringes is very high, 
and the precision on the angular diameter is reduced com- 
pared to the 102.5 m baseline UT1-UT3. 



Fig. 3. Squared visibilities obtained on £Cai on JD = 
2452742.712 (thin line) and 2452763.555 (thick line), re- 
spectively at pulsation phases 0.722 and 0.308. The two 
UD visibility models correspond to 6*tjd = 2.801 and 
3.075 mas, and take the bandwidth smearing effect into 
account. The first minimum of the visibility function (that 
never goes down to zero) occurs for baselines of approx- 
imately 199 and 181m, for an effective wavelength of 
2.18 /j,m. 



6. Angular diameters 

The object of this section is to derive the angular diam- 
eters of the Cepheids as a function of their pulsational 
phase. We discuss the different types of models that can 
be used to compute the angular diameter from the squared 
visibility measurements. 



5.6. £ Car 

As for /3Dor, the observation of £Ca,i (S = —62 deg) is 
made particularly difficult by the vignetting of the TS 
beams. Thanks to its brightness (K « 1) the SNRs are 
15-20 on the interferometric channels, and 10-15 on the 
photometric signals, using the TS and a fringe frequency 
of 242 Hz. One observation was obtained on the E0-G1 
baseline (66m ground length), and 19 measurements on 
the B3-M0 baseline. £ Car is the most observed star in our 
sample, with a total of 22 226 processed interferograms. 
Its average diameter of approximately 3 mas makes it an 
ideal target for observations with baselines of 100 to 200 m. 
On the B3-M0 baseline, we achieved projected baselines of 
89.7 to 135.0 m, corresponding to V 2 values of 8 to 42%. 
This range is ideal to constrain the visibility model and 
derive precise values of the angular diameter. 

Fig. 01 shows the squared visibility points obtained at 
two phases on £Cai. The change in angular diameter is 
clearly visible. Thanks to the variation of the projected 
baseline on sky, we have sampled a segment of the visibility 
curve. 



6.1. Uniform disk angular diameters 

This very simple, rather unphysical model is commonly 
used for interferometric studies as it is independent of any 
stellar atmosphere model. The relationship between the 
visibility V and the uniform disk angular diameter (UD) 



V(B,6 



UDJ 



2Ji(x) 



(1) 



where x = ttB #tjd /A is the spatial frequency. This func- 
tion can be inverted numerically to retrieve the uniform 
disk angular diameter #tjd- 

While the true stellar light distributions depart signifi- 
cantly from the UD model, the UD angular diameters #ud 
given in Tables to have the advantage that they can 
easily be converted to LD values using any stellar atmo- 
sphere model. This is achieved by computing a conversion 
factor #ld/$ud from the chosen intensity profile (see e.g. 
Davis et al. lSUUUI for details). 



P. Kervella et al.: VINCI/VLTI interferometric observations of Cepheids 



1.0 

0.9 
_ 0-8 

| 0.7 
E 

- 0.6 
0.5 
0.4 




1.0 0.8 0.6 0.4 0.2 0.0 

mu 



Fig. 4. Average intensity profiles computed from the four- 
parameter approximations of Claret ( 2000 ) for X Sgr (thin 
line) and t Car (thick line) , using the parameters listed in 
Tabled 

6.2. Static atmosphere intensity profile 

The visibility curve shape before the first minimum is al- 
most impossible to distinguish between a uniform disk 
(UD) and limb darkened (LD) model. Therefore, it is nec- 
essary to use a model of the stellar disk limb darken- 
ing to deduce the photospheric angular size of the star, 
from the observed visibility values. The intensity pro- 
files that we chose were computed by Claret l)2()()()|l . based 
on model atmospheres by Kurucz (1992). They consist of 
four-parameter approximations to the function 7(/i)//(l), 
where fi = cos 8 is the cosine of the azimuth of a surface 
element of the star. They are accurate approximations of 
the numerical results from the ATLAS modeling code. The 
analytical expression of these approximations is given by: 

4 

J( M )//(l) = l-5> fc (l-/z4) (2) 

k=l 

The afc coefficients are tabulated by Claret 120001) for a 
wide range of stellar parameters (T c s, logg,...) and pho- 
tometric bands (U to K). The values for each Cepheid 
are given in Tabled for the K band, and the intensity 
profiles of X Sgr and t Car are shown on Fig.^J 

The limb darkening is directly measurable by in- 
tcrfcromctry around the first minimum of the visibility 
function, as demonstrated by several authors on giant 
stars (Quirrenbach et al. 19911 Wittkowski et al. IMTT|) . 
Unfortunately, even for I Car observed in the K band, this 
requires a baseline of more than 180 m that was not avail- 
able for the measurements reported here. It is intended in 
the near future to measure directly the LD of a few nearby 
Cepheids, using the shorter wavelength bands of AMBER 
(Petrov et al. 2000 ) and the longest baselines of the VLTI 
(up to 202 m). 

6.3. Changes of limb darkening with phase 

As shown by Marengo et al. ( 2002 ) , the atmosphere of the 
Cepheids departs from that of a non-variable giant with 
identical T G ff and log g, due in particular to the presence of 
energetic shock waves at certain phases of the pulsation. 

However, this effect is enhanced at visible wavelengths 
compared to the infrared, and appears to be negligible in 
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the case of the VINCI observations. Marengo et al. ( 2003) 
have derived in the H band a relative variation of the limb 
darkening coefficient k = #tjd/#ld of only 0.2%. This is 
below the precision of our measurements and is neglected 
in the rest of this paper. Furthermore, the VINCI/VLTI 
measurement wavelength being longer (2.18/im) than the 
H band, the LD correction is even smaller, as is its ex- 
pected variation. 

From the results of Marengo et al. ( 20(321) ^ appears 
clearly that the interferometers operating at infrared 
wavelengths are ideally suited for Cepheid measurements 
that aim at calibrating the P-R and P-L relations. On the 
other hand, as pointed out by these authors, the visible 
wavelength interferometers should be favored to study the 
dynamical evolution of the atmosphere (including the limb 
darkening) during the pulsation. The geometrical determi- 
nation of the pulsation parallax is almost independant of 
the adopted atmosphere model in the K band, while this 
is not the case at shorter wavelengths. 

6.4. Visibility model and limb darkened angular 
diameters 

The VINCI instrument bandpass corresponds to the K 
band filter, transparent between A = 2.0 and 2.4 /im. An 
important effect of this relatively large spectral bandwidth 
is that several spatial frequencies are simultaneously ob- 
served by the interferometer. This effect is known as band- 
width smearing (Kervella et al.[2003b). 

To account for the bandwidth smearing, the model vis- 
ibility is computed for regularly spaced wavenumber spec- 
tral bins over the K band, and then integrated to obtain 
the model visibility. In this paper, we assume that the limb 
darkening law does not change over the K band. This is 
reasonable for a hot and compact stellar atmosphere, but 
is also coherent with the range of visibilities measured on 
the Cepheids of our sample. If necessary, this computa- 
tion can easily be extended to a wavenumber dependant 
I(/j,,a) intensity profile. Following Davis et al. {2000), us- 
ing a Hankel integral, we can derive the visibility law 
V(B, 6*ld, cr) from the intensity profile: 

V= A~J / (^) J o(^ <j0 ldv / 1-M 2 )m^ (3) 
where a is the wavenumber: 

a = 1/A (4) 
and A is a normalization factor: 

A= f I{ji)ndp (5) 
Jo 

The integral of the binned squared visibilities is com- 
puted numerically over the K band and gives the model 
V 2 for the projected baseline B and the angular diameter 
$ld through the relation: 

V 2 (9 LD ,B)= [ [V(B,9 hu ,(T)T(a)} 2 da (6) 

J K 
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where T(tr) is the normalized instrumental transmission 
defined so that 

/ T{a)dcr=l (7) 

JK 

We computed a model of T(o~) by taking into account 
the instrumental transmission of VINCI and the VLTI. 
It was first estimated by considering all known factors 
(filter, fibers, atmospheric transmission,...) and then cali- 
brated on sky based on several observations of bright stars 
with the 8 meter UTs (see Kervella et al. I2003bl for more 
details). This gives, for our sample of Cepheids, a mea- 
surement wavelength of 2.179±0.003 /im. The variation of 
effective temperature between the stars of our sample and 
over the pulsation does not change this value by more than 
±0.001 /im. The uncertainty on the effective wavelength of 
the measurement translates to a 0.15 % uncertainty on the 
measured angular diameters. Considering the level of the 
other sources of error (statistical and systematic), the ef- 
fect on our angular diameter results is negligible. 

The V 2 (9ly),B) model is adjusted numerically to the 
observed (B, V 2 ) data using a classical \ 2 minimization 
process to derive #ld- A single angular diameter is de- 
rived per observation session, the fit being done directly 
on the set of V 2 values obtained during the session. The 
systematic and statistical errors are considered separately 
in the fitting procedure, to estimate the contribution of 
the uncertainty of the calibrator diameter on the final er- 
ror bar. 

Each observation session was generally executed in less 
than 3 hours, a short time compared to the pulsation pe- 
riods of the Cepheids of our sample. Therefore, we do not 
expect any phase induced smearing from this averaging. 

6.5. Measured angular diameters 

The derived angular diameters are given in Tables to El 
for the seven Cepheids of our sample. Two error bars are 
given for each angular diameter value: 

— one statistical uncertainty, computed from the disper- 
sion of the V 2 values obtained during the observation, 

— one systematic uncertainty defined by the error bars 
on the calibrator stars a priori angular sizes. 

While the statistical error can be diminished by repeatedly 
observing the target, the systematic error is not reduced 
by averaging measurements obtained using the same cali- 
brator. 

The reference epochs To and periods P for each 
Cepheid are given in Tabled N is the number of batches 
(500 interferograms) recorded during the corresponding 
observing session. For each angular diameter, the statisti- 
cal and systematic calibration errors are given separately, 
except for the FLUOR/IOTA measurements of ( Gem, for 
which the systematic calibration error is negligible com- 
pared to the statistical uncertainty. 



7. Linear diameter curves 

For each star we used radial velocity data found in the lit- 
erature. Specifically, we collected data from Bersier (2002) 
for TyAql, ^Car, and /3Dor; from Bersier et al. (1994) for 
CGcm; from Babel et al(HHI> for WSgr. All these data 
have been obtained with the CORAVEL radial velocity 
spectrograph (Baranne, Mayor & Poncet ]T979f> . We also 
obtained data from Evans & Lyons (1986) for YOph and 
from Wilson et al. (fTTTHT^) for X Sgr. 

In theory, the linear diameter variation could be de- 
termined by direct integration of pulsational velocities 
(within the assumption that the r = 1 photosphere is 
comoving with the atmosphere of the Cepheid during its 
pulsation). However these velocities are deduced from the 
measured radial velocities by the use of a projection factor 
p. The Cepheid's radii determined from the BW method 
depend directly from a good knowledge of p. Sabbey et 
al. (P.995) and Krockenberger et al. I|1997|l have studied 
in detail the way to determine the p-factor. We used a 
constant projection factor p — 1.36 in order to trans- 
form the radial velocities into pulsation velocities. Burki 
et al. (|1982(l have shown that this value is appropriate for 
the radial velocity measurements that we used. 

8. Cepheids parameters 

8.1. Angular diameter model fitting and distance 
measurement 

From our angular diameter measurements, we can derive 
both the average linear diameter and the distance to the 
Cepheids. This is done by applying a classical \ 2 min- 
imization algorithm between our angular diameter mea- 
surements and a model of the star pulsation. The mini- 
mized quantity with respect to the chosen model is 

^2 _ V"^ (^LDobserv(^i) — #LD model {4>i ) ) 2 ^ 
^obscrv 

where cf>i is the phase of measurement i. The expression 
of #ld model {4>i) is defined using the following parameters: 

— the average LD angular diameter #ld (in mas) , 

— the linear diameter variation AD(<pi) (in -D©), 

— the distance d to the star (in pc). 

The resulting expression is therefore: 

0LDmodcl(&) = #LD~ + 9.305 W (9) 

As AD ((f>i) is known from the integration of the radial 
velocity curve (Sect.Q), the only variable parameters are 
the average LD angular diameter #ld and the distance 
d. From there, three methods can be used to derive the 
distance d, depending on the level of completeness and 
precision of the angular diameter measurements: 

— Constant diameter fit (order 0): The average lin- 
ear diameter D of the star is supposed known a pri- 
ori from previously published BW measurements or 
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JD Stations Baseline Phase #ud (mas) 6<ld (mas) V x?cd Calibrators 

(m) ± stat. ± syst. ± stat. ± syst. 



2452741.903 


B3- 


-M0 


138.366 


0.560 


1.458 ± 


0.048 


± 


0.032 


1.495 


± 


0.049 


± 


0.033 


2 


0.66 


X 


Sco 


2452742.885 


B3- 


■MO 


137.432 


0.700 


1.511 ± 


0.058 


± 


0.034 


1.549 


± 


0.059 


± 


0.035 


3 


0.52 


X 


Sco 


2452743.897 


B3- 


•MO 


137.903 


0.844 


1.415 ± 


0.055 


± 


0.034 


1.451 


± 


0.057 


± 


0.035 


3 


0.08 


X 


Sco 


2452744.868 


B3- 


■MO 


139.657 


0.983 


1.460 ± 


0.051 


± 


0.029 


1.497 


± 


0.052 


± 


0.030 


2 


0.09 


X 


Sco 


2452747.848 


B3- 


■MO 


139.530 


0.408 


1.499 ± 


0.213 


± 


0.038 


1.537 


± 


0.219 


± 


0.039 


1 




X 


Sco 


2452749.832 


B3- 


■MO 


139.084 


0.691 


1.429 ± 


0.099 


± 


0.034 


1.465 


± 


0.101 


± 


0.034 


2 


0.35 


X 


Sco 


2452766.811 


B3- 


-MO 


138.853 


0.112 


1.393 ± 


0.070 


± 


0.036 


1.428 


± 


0.071 


± 


0.037 


4 


0.09 


X 


Sco 


2452768.877 


B3- 


-MO 


128.228 


0.406 


1.413 ± 


0.016 


± 


0.028 


1.449 


± 


0.016 


± 


0.029 


6 


0.62 


X 


Sco 



Table 4. Angular diameter measurements of rj Aql. 



JD 




Stations 


Baseline 
(m) 


Phase 




6>ud (mas) 
± stat. ± syst 






6>ld (mas) 
± stat. ± syst 




N 


Xred 


Calibrators 


2452524. 


564 


E0-G1 


60.664 


0.741 


1. 


.746 ± 0.100 ±0 


.074 


1. 


.792 ± 0.103 ±0 


.076 


3 





.08 


70 Aql 




2452557. 


,546 


B3-M0 


137.625 


0.336 


1. 


.877 ± 0.098 ± 


,037 


1. 


.931 ±0.101 ±0 


,038 


1 






elnd 




2452559. 


535 


B3-M0 


138.353 


0.614 


1. 


.806 ± 0.037 ± 


,027 


1. 


.857 ± 0.038 ± 


,027 


1 






7Aqr, 


elnd 


2452564, 


532 


B3-M0 


136.839 


0.310 


1 


.809 ± 0.043 ± 


,031 


1. 


.860 ± 0.045 ± 


,032 


3 





,42 


7Aqr, 


elnd 


2452565. 


516 


B3-M0 


138.495 


0.447 


1 


.871 ±0.017 ±0. 


,027 


1. 


.924 ± 0.017 ±0 


,028 


3 





.13 


7Aqr 




2452566. 


519 


B3-M0 


137.845 


0.587 


1. 


.861 ±0.023 ±0, 


,026 


1. 


.914 ± 0.024 ±0 


,026 


5 





.23 


7Aqr 




2452567. 


523 


B3-M0 


137.011 


0.727 


1 


.802 ± 0.027 ± 


,030 


1. 


.853 ± 0.028 ±0 


,030 


2 





.62 


7Aqr 




2452573. 


511 


B3-M0 


136.303 


0.561 


1 


.884 ± 0.053 ± 0. 


,022 


1. 


.938 ± 0.054 ± 


,022 


1 






A Cm, 


HR8685 


2452769, 


937 


B3-M0 


139.632 


0.931 


1. 


.647 ± 0.026 ± 


,018 


1. 


.693 ± 0.026 ± 


,018 


3 





.06 


XSco 




2452770, 


922 


B3-M0 


139.400 


0.068 


1. 


.791 ± 0.041 ± 


,027 


1 


.842 ± 0.042 ± 


,028 


3 





,15 


XSco 




2452772, 


899 


B3-M0 


138.188 


0.343 


1. 


.880 ± 0.044 ± 


,026 


1 


.934 ± 0.046 ± 


,027 


3 





.16 


7Aqr 





Table 5. Angular diameter measurements of WSgr. 



JD Stations Baseline Phase #ud (mas) #ld (mas) N x?cd Calibrators 

(m) ± stat. ± syst. ± stat. ± syst. 



2452743.837 


B3- 


-M0 


137.574 


0.571 


1 


.408 ± 


0.096 


± 





.038 


1. 


.447 


± 


0.099 


± 





,039 


1 




X 


Sco 


2452744.915 


B3- 


-M0 


137.166 


0.713 


1. 


.292 ± 


0.088 


± 





,034 


1. 


.327 


± 


0.090 


± 





,035 


2 


0.04 


X 


Sco 


2452749.868 


B3- 


-M0 


139.632 


0.365 


1. 


.262 ± 


0.141 


± 





,040 


1. 


.297 


± 


0.145 


± 





,041 


1 




X 


Sco 


2452751.866 


B3- 


-M0 


139.538 


0.628 


1. 


.320 ± 


0.174 


± 





,041 


1. 


.357 


± 


0.179 


± 





,042 


1 




X 


Sco 


2452763.888 


B3- 


-M0 


131.830 


0.211 


1. 


.284 ± 


0.019 


± 





,029 


1. 


.319 


± 


0.020 


± 





,030 


4 


0.73 


X 


Sco 


2452764.856 


B3- 


-M0 


135.926 


0.339 


1. 


.355 ± 


0.021 


± 





,021 


1. 


.393 


± 


0.021 


± 





,022 


4 


0.76 


X 


Sco 


2452765.880 


B3- 


-M0 


132.679 


0.473 


1 


.313 ± 


0.022 


± 





,025 


1. 


.349 


± 


0.023 


± 





,026 


4 


1.43 


X 


Sco 


2452767.867 


B3- 


-M0 


132.637 


0.735 


1. 


.208 ± 


0.073 


± 





,039 


1. 


.241 


± 


0.075 


± 





,040 


3 


0.01 


X 


Sco 


2452769.914 


B3- 


-M0 


120.648 


0.005 


1. 


.240 ± 


0.055 


± 





,034 


1. 


.274 


± 


0.056 


± 





,035 


2 


0.33 


X 


Sco 



Table 6. Angular diameter measurements of /3Dor. 



JD Stations Baseline Phase #ud (mas) 6<ld (mas) N x?cd Calibrators 

(m) ± stat. ± syst. ± stat. ± syst. 



2452215.795 


U1-U3 


89.058 


0.161 


1 


.842 


± 


0.036 


± 


0.074 


1.896 


± 


0.036 


± 


0.074 


3 


0.03 


X Phe, 7 2 Vol 


2452216.785 


U1-U3 


89.651 


0.261 


1. 


.954 


± 


0.026 


± 


0.040 


2.011 


± 


0.026 


± 


0.040 


7 


0.10 


7 2 Vol 


2452247.761 


U1-U3 


83.409 


0.408 


1. 


.921 


± 


0.045 


± 


0.039 


1.977 


± 


0.045 


± 


0.039 


5 


0.40 


eRet 


2452308.645 


U1-U3 


75.902 


0.594 


1. 


.844 


± 


0.027 


± 


0.071 


1.897 


± 


0.027 


± 


0.071 


5 


1.01 


HD 63697 


2452567.827 


B3-M0 


134.203 


0.927 


1. 


.793 


± 


0.039 


± 


0.049 


1.848 


± 


0.039 


± 


0.049 


1 




HR2549 


2452744.564 


B3-M0 


89.028 


0.884 


1. 


.730 


± 


0.064 


± 


0.032 


1.780 


± 


0.064 


± 


0.032 


2 


0.09 


HR3046, 4831 


2452749.514 


B3-M0 


98.176 


0.387 


1. 


.921 


± 


0.106 


± 


0.029 


1.978 


± 


0.106 


± 


0.029 


3 


0.11 


HR3046 


2452750.511 


B3-M0 


98.189 


0.488 


1 


.864 


± 


0.065 


± 


0.039 


1.919 


± 


0.065 


± 


0.039 


2 


0.24 


HR3046 


2452751.519 


B3-M0 


95.579 


0.591 


1 


.954 


± 


0.169 


± 


0.030 


2.012 


± 


0.169 


± 


0.030 


3 


0.03 


HR3046 
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Table 7. VINCI/VLTI and FLUOR/IOTA angular diameter measurements of £ Gem. No systematic calibration error 
is given for FLUOR/IOTA values (negligible compared to the statistical uncertainty). The baseline is given for the 
VINCI/VLTI observations (in m), while the spatial frequency (in italic) is listed for the measurements obtained with 
FLUOR, expressed in cycles/arcsec. 



JD 


Stations 


B, SF 


Phase 


6» UD (mas) 


<9 LD (mas) 


N 


Xrcd 


Calibrators 










± stat. ± syst. 


± stat. ± syst. 








2452214.879 


U1-U3 


82.423 


0.408 


1.677 ±0.030 ±0.051 


1.725 ± 0.031 ± 0.052 


8 


0.25 


39 Eri 


2452216.836 


U1-U3 


72.837 


0.600 


1.712 ±0.057 ±0.067 


1.760 ±0.058 ±0.069 


6 


0.28 


39 Eri, 7 2 Vol 


2451527.972 


IOTA-38m 


81870 


0.739 


1.606 ± 0.334 


1.651 ±0.343 


1 




HD 49968 


2451601.828 


IOTA-38m 


83.917 


0.014 


1.709 ± 0.086 


1.795 ± 0.088 


3 


0.02 


HD 49968 


2451259.779 


IOTA-38m 


83.760 


0.318 


2.040 ± 0.291 


2.144 ± 0.299 


1 




HD 49968 


2451262.740 


IOTA-38m 


84.015 


0.610 


1.692 ± 0.273 


1.767 ±0.281 


2 


0.13 


HD 49968 


2451595.863 


IOTA-38m 


83. 790 


0.427 


1.391 ± 0.284 


1.306 ± 0.292 


2 


1.72 


HD 49968 


2451602.764 


IOTA-38m 


85.010 


0.107 


1.867 ±0.216 


1.962 ± 0.222 


2 


0.02 


HD 49968 



Table 8. Angular diameter measurements of YOph. 



JD Stations Baseline Phase #ud (mas) 6>ld (mas) N Xrcd Calibrators 

(m) ± stat. ± syst. ± stat. ± syst. 



2452742.906 


B3-M0 


139.569 


0.601 


1.427 


± 


0.115 


± 


0.034 


1.472 


± 


0.119 


± 


0.035 


2 


0.10 


X 


Sco 


2452750.884 


B3-M0 


139.057 


0.067 


1.380 


± 


0.100 


± 


0.034 


1.423 


± 


0.103 


± 


0.035 


2 


0.41 


X 


Sco 


2452772.831 


B3-M0 


139.657 


0.349 


1.443 


± 


0.051 


± 


0.025 


1.488 


± 


0.053 


± 


0.026 


3 


0.22 


X 


Sco 


2452782.186 


B3-M0 


129.518 


0.168 


1.402 


± 


0.027 


± 


0.037 


1.445 


± 


0.028 


± 


0.038 


4 


0.30 


X 


Sco 



Table 9. Angular diameter measurements of ^Car. 



JD Stations Baseline Phase #ud (mas) #ld (mas) N Xrcd Calibrators 

(m) ± stat. ± syst. ± stat. ± syst. HR 



2452453. 


.498 


E0-G1 


61.069 





.587 


2 


.958 


± 





.039 


± 





.102 


3. 


.054 


± 





.041 


± 





.105 


4 





.01 


4050 




2452739. 


.564 


B3-M0 


130.468 





.634 


2 


.786 


± 





.073 


± 





.042 


2. 


891 


± 





.076 


± 





.043 


2 





.03 


4526 




2452740. 


.569 


B3-M0 


128.821 





.662 


2 


.879 


± 





.017 


± 





.042 


2. 


.989 


± 





.018 


± 





.044 


7 





.77 


4526 




2452741. 


.717 


B3-M0 


96.477 





.694 


2 


.893 


± 





.025 


± 





.028 


2. 


993 


± 





.026 


± 





.029 


5 





.28 


4526 




2452742. 


.712 


B3-M0 


99.848 





.722 


2 


.801 


± 





.034 


± 





.042 


2. 


.899 


± 





.035 


± 





.043 


5 





.09 


4526 




2452743. 


.698 


B3-M0 


99.755 





.750 


2 


.667 


± 





.071 


± 





.015 


2. 


.758 


± 





.074 


± 


0. 


.016 


2 


0. 


.08 


4831 




2452744. 


634 


B3-M0 


114.981 





.776 


2 


.698 


± 





.031 


± 





.012 


2. 


.794 


± 





.032 


± 





.013 


6 





.73 


4831 




2452745. 


.629 


B3-M0 


115.791 





.804 


2 


.584 


± 





.094 


± 





.017 


2. 


.675 


± 





.097 


± 





.017 


2 





.01 


3046, 


4546, 4831 


2452746. 


.620 


B3-M0 


116.828 





.832 


2 


.679 


± 





.023 


± 





.039 


2. 


.775 


± 





.023 


± 





.040 


5 





.65 


3046, 


4546 


2452747. 


.599 


B3-M0 


120.812 





.860 


2 


.606 


± 





.122 


± 





.025 


2. 


.699 


± 





.127 


± 





.026 


3 


0. 


.70 


4546, 


4831 


2452749. 


.576 


B3-M0 


124.046 





.915 


2 


.553 


± 





.075 


± 





.011 


2. 


.645 


± 





.077 


± 


0. 


.012 


4 


1. 


.18 


4546 




2452751. 


.579 


B3-M0 


122.555 





.971 


2 


.657 


± 





.027 


± 





.017 


2. 


.753 


± 





.028 


± 





.017 


4 


1. 


.16 


3046, 


4831 


2452755. 


.617 


B3-M0 


112.185 





.085 


2 


.867 


± 





.109 


± 





.013 


2. 


.970 


± 





.113 


± 


0. 


.013 


1 






4831 




2452763. 


.555 


B3-M0 


120.632 





.308 


3 


.077 


± 





.008 


± 





.031 


3. 


194 


± 





.009 


± 





.033 


6 


1 


.02 


4546 




2452765. 


.555 


B3-M0 


119.629 





.365 


3 


.094 


± 





.011 


± 





.031 


3. 


.212 


± 





.011 


± 


0. 


.033 


6 


1 


.19 


4546 




2452766. 


.550 


B3-M0 


120.005 





.393 


3 


.092 


± 





.011 


± 





.032 


3. 


210 


± 





.011 


± 


0. 


.033 


7 


0. 


.99 


4546 




2452768. 


.566 


B3-M0 


115.135 





.450 


3 


.075 


± 





.010 


± 





.034 


3. 


.188 


± 





.011 


± 





.035 


7 


0. 


.46 


4546 




2452769. 


.575 


B3-M0 


113.082 





.478 


3 


.075 


± 





.018 


± 





.011 


3. 


189 


± 





.018 


± 


0. 


.012 


3 





.03 


3046, 


4831 


2452770. 


535 


B3-M0 


121.152 





.505 


3 


.044 


± 





.019 


± 





.009 


3. 


160 


± 





.020 


± 





.009 


2 





.20 


3046, 


4831 


2452771. 


.528 


B3-M0 


122.014 





.533 


3 


.021 


± 





.017 


± 





.010 


3. 


136 


± 





.017 


± 





.010 


3 





.88 


4831 





P-R relations (See Sect. 18. 2[) . We assume here that 
AD(<p) — 0. The only remaining variable to fit is the 
distance d. This is the most basic method, and is use- 
ful as a reference to assess the level of detection of the 
pulsational diameter variation with the other methods. 
— Variable diameter (order 1): We still consider that 
the average linear diameter D of the star is known a 
priori, but we include in our model the radius varia- 



tion derived from the integration of the radial veloc- 
ity curve. This method is well suited when the intrin- 
sic accuracy of the angular diameter measurements is 
too low to measure precisely the pulsation amplitude 
(C Gem, X Sgr and Y Oph) . The distance d is the only 
free parameter for the fit. 
— Complete fit (order 2): The average LD angular 
diameter 6*ld and the distance d are both considered 
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as variables and adjusted simultaneously to the an- 
gular diameter measurements. In the fitting process, 
the radius curve is matched to the observed pulsation 
amplitude. Apart from direct trigonometric parallax, 
this implementation of the BW method is the most di- 
rect way of measuring the distance and diameter of a 
Cepheid. It requires a high precision angular diameter 
curve and a good phase coverage. It can be applied 
directly to our ?yAql, WSgr, /3Dor and £Car measure- 
ments. 



8.2. Published linear diameter values 

In this section, we survey the existing linear diameter de- 
terminations for the Cepheids of our sample, in order to 
apply the order and 1 methods to our observations. 

A large number of BW studies have been published, 
using both visible and infrared wavelength observations. 
For £Gem and 7?Aql, the pulsation has been resolved us- 
ing the Palomar Testbed Interferometer (Lane et al.|2000 
2002), thererefore giving a direct estimate of the diam- 
eter and distance of these stars. TablclTUl gives a list of 
the existing diameter estimates for the Cepheids of our 
sample from the application of the classical BW method 
("Baade-Wesselink" section of the table). 

From the many different P-R relations available, we 
chose the Gieren et al. Ijl998[l version, as it is based on in- 
frared colors for the determination of the temperature of 
the stars. Compared to visible colors, the infrared colors 
give a much less dispersed P-R relation. Indeed, this re- 
lation has a very good intrinsic precision of the order of 5 
to 10 % for the period range of our sample. Moreover, it is 
identical to the law determined by Laney & Stobic ( 1995). 
The compatibility with the individual BW diameter esti- 
mates is also satisfactory. The linear diameters deduced 
from this P-R law are mentioned in the "Empirical P- 
R" section of TableHIH We assume these linear diameter 
values in the following. 

8.3. Angular diameter fitting results 

The results of both constant and variable diameter fits for 
the seven Cepheids of our sample are listed in TablclTTl 
to 1131 r?Aql, WSgr, /3Dor and £Car gave results for all 
fitting methods, while X Sgr, £ Gem and Y Oph were lim- 
ited to order 1 models. For XSgr, the order 1 fit is less 
adequate than the order 0, considering the quality of our 
measurements of this star. This is shown by the fact that 
the x 2 is significantly higher for the order 1 fit (1.36) than 
for the order (0.38). 

In the case of I Car, the fit of a constant diameter re- 
sults in a very high \ 2 value. This means that the average 
diameters #udo and 6>ldo should not be used for further 
analysis. The pulsation curve of this star is not sampled 
uniformly by our interferometric observations, with more 
values around the maximum diameter. This causes the 
larger diameter values to have more weight in the aver- 



1.80 



% 1.70 
E 



I 1.60 



, 1.50 



1.40 ~- 



1.30 



1.20 



0.0 



0.2 



0.4 



1.0 



Phase 



Fig. 5. Order model fit for XSgr. 



2.10 




o 1.70 - 



1.60 - 



0.0 



0.2 



0.4 




Phase 



Fig. 6. Order 2 model fit for rj Aql. The superimposed an- 
gular diameter variation curve (thin line) is derived from 
the integration of the radial velocity curve. 



age diameter computation, and this produces a significant 
positive bias. This remark does not apply to the orders 1 
and 2 fitting methods. 

As a remark, no significant phase shift is detected at 
a level of 2.5 10~ 4 (14 minutes of time) between the pre- 
dicted radius curve of t Car and the observed angular di- 
ameter curve. The values of P and To used for the fit are 
given in Table 

Fig. El to ^2 show the best models for each star, to- 
gether with the VINCI/VLTI angular diameter measure- 
ments for the seven Cepheids of our sample. Fig.H2l gives 
an enlarged view of the maximum diameter of £ Car. 
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Table 10. Published linear diameter estimates, expressed in Dq. 




XSgr 


rj Aql 


W Sgr 


(5 Dor 


C Gem 


Y Oph 


I Car 


Interferometry 
















Kervella et al. «2001bl* 
















Lane et al. f2002H 




61.8 ±7.6 






66.7 ±7.2 






Nordgren et al. ( 2000)* 










60±? 4 






Baade-Wesselink 
















Bersier et al. (I1997H 






56.0 ±2.9 




89.5 ± 13.3 






Fouque et al. (ISTMt 




48.1 ± 1.1 










201.7 ±3.0 


Krockenberger et al. 1199711 






56.8 ±2.3 




69.111J 






Laney & btobie ( 199o ) 








63.0 ± l.o 




92.2 ± 6.1 


i on i _l_ A c 
loO.l ± 4.0 


Moffett & Barnes ltTWIl a 
MotTett & Barnes ltTWt b 


47.8 ± 4.5 


52.8 ±3.8 


60.8 ± 7.6 




62.6 ± 11.5 






49.6 ±4.6 


54.8 ±3.9 


63.1 ±7.8 




64.9 ± 11.9 






Sabbey et al. tTM& c 


42.2 ±4.1 


62.7 ±3.1 






61.8 ±3.5 






Sabbev et al. Ji99Sl d 


66.6 ± 4.9 


65.8 ± 3.2 






64.4 ± 3.6 






Sachkov et al. itTTOl 










74 ± 10 






Taylor et al. ftWl 














179.2 ± 10.4 


Taylor & Booth RMil 








67.8 ±0.7 








Turner & Burke 11301721 




52.6 ±8.9 




53.8 ±1.9 








Sasselov & Lester H1990t 


67 ±6 


62 ±6 












Mean B-W (overall a) 


52.5 (11.4) 


59.9 (5.7) 


57.0 (3.4) 


65.8 (7.2) 


65.3 (9.8) 


92.2 (-) 


180 (-) 


Empirical P-R 
















Gieren et al. (11998 ll 


51.2 ±2.6 


52.1 ± 2.7 


54.4 ±2.9 


66.0 ±3.9 


67.6 ±4.0 


100.1 ±7.3 


173.1 ± 15.8 



* £Gem values were derived from Kervella et al. (2001b) and Nordgren et al. (2000) using the HlPPARCOS parallaxes. 77 Aql 

was taken from Nordgren et al. ( 2000 ) 
a Assuming a constant p-factor. 
b Assuming a variable p-factor. 
c Bisector method. 
d Parabolic fit method. 



Table 11. Order 0. Cepheid average angular diame- 
ters and distances derived from the VINCI interferomet- 
ric measurements, assuming a constant diameter model 
(AD = 0), The average diameter D is taken from Gieren 
et al. I)1998|l . Two error bars are given in brackets for the 
angular diameter: the statistical dispersion and the cal- 
ibration systematics. The uncertainty mentioned for the 
distance d is the quadratic sum of the statistical, calibra- 
tion and P-R a priori diameter errors, the last two being 
systematic in nature. The three types of errors are also 
reported separately in brackets. The results for I Car are 
mentioned only for completeness, but are not meant to be 
used for further analysis, as our observations are inconsis- 
tent with a constant diameter model. 



Star 


#ldo (mas) 


d (pc) 


xl 


XSgr 


1.471 ± 0.033 [0 .oi 3 0.031] 


324 ± 18[ 3 7 17] 


0.38 


77 Aql 


1.856 ± 0.028[o.oo9 0.026] 


261 ± 14[i 4 14] 


3.98 


WSgr 


1.348 ± 0029 [0 .on 0.027] 


376 ± 22 [3 8 21] 


0.90 


/3Dor 


1.926 ± 0.024 [0 . 014 o. 020] 


319 ± 20 [3 2 19] 


1.31 


CGem 


1.747 ±0.061 [0 .025 0.056] 


360 ± 25 [5 12 22] 


0.51 


YOph 


1.459 ±0.040 [0 . 23 0.033] 


638 ± 50[io 14 47] 


0.16 


(£Car) 


3.071 ± 0.012[o.oo4 0.011] 


524 ± 49[i 2 49] 


23.2 



1.60 



SS 1.50 - 

E. 

CD 




I 1.10 - 



1.00 -I 1 1 1 1 

0.0 0.2 0.4 0.6 0.8 1.0 

Phase 

Fig. 7. Order 2 model fit for W Sgr. 
9. Discussion 

9.1. Limb darkening of r] Aql and ( Gem 

From the NPOI (Armstrong et al. I21M1 Nordgren et 
al. l2UUU|l . PTI (Lane et al. lSUUSjl and VINCI/VLTI mea- 
surements, we know the average UD angular diameters 
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2.20 




1.70 



1.60 - 



t 1-30 - 



1.20 




Fig. 8. Order 2 model fit for /3Dor. 



Fig. 10. Order 1 model fit for YOph. 
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1 3.00 



I 2.80 
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0.0 
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1.0 



Fig. 9. Order model fit for £ Gem. The crosses represent 
the FLUOR/IOTA data, and the two points are UT1-UT3 
observations with VINCI. 



of r\ Aql and £ Gem at several effective wavelengths with 
high precision. Table ITU gives the angular diameter values 
and the corresponding wavelengths. Claret's (2000) linear 
limb darkening parameters u were used to compute the 
expected conversion factors p — #ld /S'ud ■ To read the u 
table, we have considered the closest parameters to the 
average values for r\ Aql and £ Gem in Table ^ an d we 
computed p using the formula from Hanbury Brown el 
al. (JTHZU): 



' l-u/3 
1 - 7u/15 



(10) 



For the NPOI observation (A e s ~ 0.73 pm), we have cho- 
sen an intermediate value of u between the R and / bands. 

We note that the value of #ld for r\ Aql that we derive 
for the NPOI observation, #ld = 1.73 ± 0.04 mas, is not 
identical to the LD angular diameter originally given by 



Fig. 11. Order 2 model fit for £Car. 

Armstrong et al. lj2UUT|) . 6» LD = 1.69 ± 0.04 mas. There is 
a la difference, that may be due to the different source 
of limb darkening coefficient that these authors used for 
their modeling (Van Hamme 1993). 

The resulting 6*ld values for the three observations are 
compatible at the 2 a level, but there is a slight trend that 
points towards an underestimation of the limb darkening 
effect at shorter wavelengths, or alternatively its overesti- 
mation at longer wavelengths. Considering that the limb 
darkening is already small in the infrared, the first hy- 
pothesis seems more plausible. Marengo et al. f2003 ) 
have shown that the Cepheids limb darkening can be sig- 
nificantly different from stable giant stars, particularly at 
visible wavelengths. This could explain the observed differ- 
ence between the 0.73 pm and K band diameters of rj Aql 
and £Gem, the latter being probably closer to the true 
LD diameters, thanks to the lower limb darkening in the 
infrared. 

In the case of rj Aql, another explanation could be that 
the measurement at visible wavelengths is biased by the 
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Fig. 12. Detail of Fig. ^2 showing the angular diameter 
curve of I Car around the maximum diameter. 

Table 12. Order 1. Cepheid angular diameters and dis- 
tances, assuming the average diameter D of Gieren et al. 
(1998$. The diameter variation curve AD(<j)) is integrated 
from the radial velocity curve. Only the distance is ajusted 
by the fitting procedure. The error bars on d are given as 
in Table CH 



Star 


#ldi (mas) 


di (pc) 


2 
Xi 


XSgr 


1.461 ± 0.033[ .oi3 0.03i] 


326 ± 18(3 7 17] 


1.36 


77 Aql 


1.839 ± 0.028[o.oo9 0.026] 


264 ± % 4 i4] 


0.40 


WSgr 


1.312 ±0029,0.011 0.027] 


386 ± 22(3 8 21] 


0.42 


/3Dor 


1.884 ± 0.024 [0 .oi4 0.020] 


326 ± 20(3 2 19] 


0.23 


(Gem 


1.718 ± 0.061(0.025 0.056] 


366 ± 25(5 12 22] 


0.88 


YOph 


1.437 ± 0.040[o.o23 0.033] 


648 ±51 [10 15 47] 


0.03 


TZCar 


2.977 ± 0.012[o.oo4 0.011] 


542 ± 49[i 2 49] 


0.71 



blue companion of 77 Aql. However, it is 4.6 magnitudes 
fainter than the Cepheid in the V band (Bohm-Vitense & 
Proffitt 1985, see also Sect. l9.2(l . and therefore should not 
contribute significantly to the visibility of the fringes. 

Table 13. Order 2. Cepheid average angular diameters 
and distances determined through the application of the 
modified BW method. The only input is the diameter vari- 
ation curve AD ((f)) derived from the integration of the ra- 
dial velocity. The distance and average angular diameter 
are ajusted simultaneously. The statistical and systematic 
errors on d are listed separately in brackets. 



Star LD2 (mas) d 2 (pc) xl 



77 Aql 


1.839 ± 0.028[o.oo9 0.026] 


276±£ 


r55 61 
38 4j 


0.43 


WSgr 


1.312 ± 0.029,0.011 0.027] 


o 7 o+216 
o,y -130 


r216 111 
ll30 7 ] 


0.48 


/3Dor 


1.891 ± 0.024,0.014 0.020] 


345t£ 8 


rl75 Si 
[80 2 ] 


0.25 


^Car 


2.988 ± 0.012[o.oo4 0.011] 


603t^ 


(24 31 
19 2j 


0.49 



Table 14. Average UD angular diameter of 77 Aql and 
£ Gem from the litterature, and the associated conver- 
sion factor p = #ld/$ud from the linear limb darken- 
ing coefficients of Claret (2000 ). References: (1) Armstrong 
et al. (EHHJ) and Nordgren et al. pltHl)) . (2) Lane et 
al. (USUI), (3) this work. 



Ref. 


A (jim) 


6» UD (mas) 


P 


9 LO (mas) 


77 Aql 










(1) 


0.73 


1.65 ±0.04 


1.048 


1.73 ±0.04 


(2) 


1.65 


1.73 ±0.07 


1.024 


1.77 ±0.07 


(3) 


2.18 


1.80 ±0.03 


1.021 


1.84 ±0.03 


CGem 










(1) 


0.73 


1.48 ±0.08 


1.051 


1.56 ±0.08 


(2) 


1.65 


1.61 ±0.03 


1.027 


1.65 ±0.03 


(3) 


2.18 


1.70 ± 0.06 


1.023 


1.75 ±0.06 



9.2. Binarity and other effects 

As demonstrated by several authors (see Szabadosp003 
for a complete database), binarity and multiplicity are 
common in the Cepheid class. Evans (jl992l has observed 
that 29 % of the Cepheids of her sample have detectable 
companions. 

Our sample of Cepheids contains four confirmed bi- 
nary Cepheids, out of a total of seven stars. As it is 
biased towards bright and nearby Cepheids, this large 
fraction is an indication that many Cepheids currently 
believed to be single could have undetected compan- 
ions. XSgr (Szabados I1989bjl . 77 Aql (Bohm-Vitense & 
Proffitt EH^and WSgr (Bohm-Vitense & Proffitt 115551 
Babel et al. 1989) are confirmed members of binary or 
multiple systems. £ Gem is a visual binary star (Proust 
et al. I1981fl . but the separated companion does not con- 
tribute to our observations. Y Oph was once suspected to 
be a binary (Pel I1978|l . but Evans (|1992|l has not con- 
firmed the companion, and has set an upper limit of AO 
on its spectral type. 

The physical parameters of the companions of 77 Aql 
and WSgr have been derived by Bohm-Vitense & 
Proffitt (1985) and Evans l|1991fl . based on ultravio- 
let spectra. The latter has derived spectral types of 
B9.8V and A0V, respectively. The orbital parameters of 
the binary W Sgr were computed by Babel et al. l)1989f) 
and Albrow & Cottrell ( 1996 ). Based on IUE spectra, 
Evans (1992) has set an upper limit of AO on the spec- 
tral type of the companion of X Sgr. 

The difference in V magnitude between these three 
Cepheids and their companions is AMy > 4.5. The AMk 
is even larger due to the blue color of these stars, AMk > 
5.7. Therefore, the effect on our visibility measurements 
is negligible, with a potential bias of AV 2 < 0.5 %. For 
example, this translates into a maximum error of ±11 pas 
on the average angular diameter of 77 Aql, (a relative er- 
ror of ±0.6 %), that is significantly smaller than our error 
bars (±1.5%). In the K band, the effect of the compan- 
ions of the other Cepheids is also negligible at the pre- 
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cision level of our measurements. However, the presence 
of companions will have to be considered for future mea- 
surements with angular diameter precisions of a few fj,as. 
In this respect, long-period Cepheids, such as ^Car, are 
more reliable, as their intrinsic brightness is larger than 
the short-period pulsators, and therefore they dominate 
their potential companions even more strongly. 

Fernie et al. (1995b) have found that the amplitude 
of the light curve of Y Oph has been decreasing for a few 
decades. A similar behavior has been observed only on 
Polaris (e.g. Evans et al. 2002). The uncertainty on our 
#ld measurements has not allowed us to detect unam- 
biguously the pulsation of this star, but it is clearly an im- 
portant target for future observations using the Auxiliary 
Telescopes (1.8 m) of the VLTI in order to estimate its 
parameters with high precision. 

Interestingly, Gieren et al. Q1993[) have studied the im- 
pact of binary Cepheids on their determination of the 
period-luminosity relation using 100 Cepheids, and they 
conclude that it is negligible. This is due to the very large 
intrinsic luminosity of the Cepheids that overshine by sev- 
eral orders of magnitude most of the other types of stars. 

10. Conclusion and perspectives 

We have reported in this paper our long-baseline inter- 
ferometric observations of seven classical Cepheids using 
the VINCI/VLTI instrument. For four stars (77 Aql, W Sgr, 
(3 Dor and I Car), we were able to apply a modified version 
of the BW method, resulting in an independent estimate 
of their distance. For all stars, we also derived their dis- 
tances from lower order fitting methods, that use an a 
priori estimate of their linear diameter from the P-R re- 
lation of Gieren et al. (lOSHIl. We would like to emphasize 
that the order 0/1 and order 2 error bars are different in 
nature, and they should be treated differently in any fur- 
ther use of these results. While the order 2 error bars can 
be treated as statistical (i.e. reduced by averaging), the 
order 0/1 methods errors are dominated by the system- 
atic uncertainty introduced by the a priori estimation of 
the linear radius. The respective contributions of the sta- 
tistical and systematic uncertainties are given separately 
in Tables ITTI and IT5| These values assume a constant value 
of the p-factor of 1.36, and can be scaled linearly for other 
values. 

We will use these distances in Paper II, together with 
previously published measurements, to calibrate the zero 
points of the Period-Radius and Period-Luminosity rela- 
tions In Paper III, we will calibrate the surface brightness- 
color relation, with a particular emphasis on the evolution 
of £Car in this diagram over its pulsation. These three 
empirical relations are of critical importance for the ex- 
tragalactic distance scale. 

The direct measurement of the limb darkening of 
nearby Cepheids by interferometry is the next step of the 
interferometric study of these stars. It will allow a refined 
modeling of the atmosphere of these stars. This obser- 
vation will be achieved soon using in particular the long 



baselines of the VLTI equipped with the AMBER instru- 
ment, and the CHARA array for the northern Cepheids. 
Another improvement of the interferometric BW methow 
will come from radial velocity measurements in the near 
infrared (see e.g. Butler & Bell 02021 . They will avoid 
any differential limb darkening between the interferomet- 
ric and radial velocity measurements, and therefore make 
the resulting distances more immune to limb darkening 
uncertainties. 
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